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Abstract 

The possibility of maintaining entanglement in a quantum system at finite, even 
high, temperatures - the so-called ‘hot entanglement’ - has obvious practical 
interest, but also requires closer theoretical scrutiny. Since quantum entangle¬ 
ment in a system evolves in time and is continuously subjected to environmental 
degradation, a nonequilibrium description by way of open quantum systems is 
called for. To identify the key issues and the contributing factors that may per¬ 
mit ‘hot entanglement’ to exist, or the lack thereof, we carry out a model study 
of two spatially-separated, coupled oscillators in a shared bath depicted by a 
finite-temperature scalar field. From the Langevin equations we derived for the 
normal modes and the entanglement measure constructed from the covariance 
matrix we examine the interplay between direct coupling, field-induced interac¬ 
tion and finite separation on the structure of late-time entanglement. We show 
that the coupling between oscillators plays a crucial role in sustaining entan¬ 
glement at intermediate temperatures and over finite separations. In contrast, 
the field-induced interaction between the oscillators which is a non-Markovian 
effect, becomes very ineffective at high temperature. We determine the critical 
temperature above which entanglement disappears to be bounded in the lead¬ 
ing order by the inverse frequency of the center-of-mass mode of the reduced 
oscillator system, a result not unexpected, which rules out hot entanglement in 
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such settings. 
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1. Introduction 

Recently Galve et al mia pointed out the possibility of keeping quantum 
entanglement alive in a system at high temperatures by driving the system 
of two oscillators with a time-dependent interaction term. This is important 
in practical terms because if entanglement in a quantum open system can be 
maintained at high temperatures, it eases the way how devices for quantum 
information processing can be conceptualized and designed. From a theoretical 
viewpoint understanding the basic mechanisms of obtaining this so-called ‘hot 
entanglement’ [3] is also of great interest. 

Before beginning the analysis, we note the word ‘hot’ conveys three layers 
of meaning in three different contexts, referring to quantum systems A) kept in 
thermal equilibrium at all times, B) in a nonequilibrium condition and evolving, 
possibly but not necessarily, toward an equilibrium state, and C) in a nonequilib¬ 
rium steady state at late times. In this study we derive the fully nonequilibrium 
dynamics of a system of two coupled quantum harmonic oscillators interact¬ 
ing with a common bath described by a bosonic field at finite temperature T. 
Thus our present work falls under Case B, which is in contrast to Case A 015 ], 
where a quantum system is assumed to be already in equilibrium and remains 
that way. We depict how entanglement of the open quantum system evolves 
in time and derive the critical temperature above which entanglement cannot 
survive. In an accompanying paper [6] we study one subcase of Case C, that of a 
quantum system in nonequilibrium steady state (NESS) at late times, using the 
framework and results obtained in 0. The system we analyze there consists of 
two coupled quantum harmonic oscillators each interacting with its own bath, 
described by a scalar field, set at two different temperatures T) > T 2 which 
together form the environment. Carrying out a fully systematic analysis of how 
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quantum entanglement in open systems under different nonequilibrium condi¬ 
tions evolves is, in our view, a necessity before any claim of “hot entanglement” 
can be asserted. 

In terms of methodology our present study makes use of the conceptual 
framework of quantum open systems 0 and the techniques of nonequilibrium 
quantum field theory [9]. It is a finite temperature generalization of our recent 
work m where the entanglement behavior at late times between two coupled 
and spatially separated oscillators interacting with a common bath modeled by 
a scalar field at zero temperature is analyzed in detail. That work in turn is a 
generalization of the paper of Lin and Hu m with coupling between the two 
oscillators added in the consideration. 


2. System Setup 


Our system is made up of two spatially separated coupled detectors, which 
are entities with internal degrees of freedom (idf) Xi, 2 - The idf of each detector 
is described by a harmonic oscillator of mass m and bare frequency uib- This 
system is placed in a common finite-temperature bath modeled by a massless 
scalar field <j> initially prepared in a thermal state at temperature /3 _1 . The 
system is allowed to interact with the bath initially at t = 0. We want to 
track down its evolution in time, derive the entanglement dynamics between 
the two detectors at late times and determine the critical temperature above 
which entanglement no longer exists. 

The action of the whole system is 

S[x,4>] = Jds y x?( s ) ~ ~y~ *i( s ) _ J ds ma xi(s)x 2 («) 

J d 4 x j(x)(j>(x) + j d 4 x ^ d^d^cj ), (1) 


L 2=l 


+ 


where the current j( x) takes the form j(x) = e J2f=i Xi(t) ^^ 3 ' 1 [ x — The 

spacetime coordinate x is understood as a shorthand notation of (t,x). The 
parameter er in the action is the coupling strength between the two idfs, while e 
is the coupling constant between each idf and the bath. We have written down 
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the action to allow for the detectors to move along an arbitrary yet prescribed 
trajectory z i(t). In this work we assume they stay at rest throughout. 

When the initial state of the idf has a Gaussian form, the reduced density 
matrix of the idf can be found exactly with the help of the influence functional 
formalism in the closed-time path integral framework. This enables us to obtain 
the full-time dynamics of the reduced system under the influence of the environ¬ 
ment for arbitrary coupling strengths, as was done in full detail in OH- Here, 
to highlight the physics behind thermal entanglement, we opt for a simpler, 
more physically transparent yet no less general way, by means of the Langevin 
equation approach, which has been shown to be totally compatible with the 
reduced-density-matrix description for linear systems. For the current configu¬ 
ration, the Langevin equations of, say, xi is given by 

mxi(t) + mul xi (*) + ma x 2 (i) 



In Eq. |2]), in addition to the restoring force —mw 2 x 1 an d the direct cou¬ 
pling max 2 (t) with the other idf, the essential (most interesting) physics is 
contained in the nonlocal interactions generated by the system’s interaction 
with its environment, and the stochastic driving force £1 which recounts both 
the quantum and thermal noises originating from the heat bath at the loca¬ 
tion of Detector 1. It obeys the Gaussian statistics with (£i(t)) = 0 and 


1 
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(£i(*)£i(0) = e 2 G ff (zi,t;zi,t'), where G H (x,x') 


{{4>{x),<t>(x')}), wi th { ,} 


denoting symmetrization, is the Hadamard function of the scalar field. In addi¬ 
tion, nonzero correlation of the bath between the locations of detector 1 and 2 
implies (£1 (t)^ 2 (t')) = e 2 G# (z 1; f; z 2 , t'). The (• • •) can represent the ensemble 
average or the quantum expectation values, depending on the context. 


The nonlocal expressions in ([ 2 ]) containing the retarded Green function 


Gr(x,x’) = i6(t — t')[(f>(x), 4>(x')] of the scalar field, with [,] denoting anti- 
symmetrization, embrace the dissipative self-force and the history-dependent 
non-Markovian interaction between the two idfs as the consequences of cou¬ 
pling between the idfs and the bath. In particular, these nonlocal expressions 
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are independent of the initial bath state. Essentially the stochastic forcing term 
and the nonlocal terms in |2]) capture the overall influences from the environ¬ 
ment. The temporal Fourier transforms of these two kernel functions Gr and 
Gr are connected via the fluctuation-dissipation relation, 

Bk — r°° cIk — 

Gif(R, k) = coth — ImGfl(R, k), where G(R,r)= / —G(R, k) e _iKT . 

2 J-oo 2 tt 

(3) 

In certain contexts it signifies a balance between the energy transfer via noise 
from, and the dissipation back to, the environment. Thus the stochastic equa¬ 
tions of motion of Xi, Xi describe a set of coupled, damped, driven oscillators 
undergoing non-Markovian dynamics. 

3. Dynamics 

The set of equations of motion for xi, X 2 in fact can be decoupled into the 
center of mass (CoM) mode x+ = (xi + X2)/2 and the relative mode X- = 
Xi - X 2 HD], 

X+ (i) + 2 7 x+ (t) - 2 7 6<yt ^ x+ (t - i) + w+ x+ (*) = —£+ (t), (4) 

i m 

x- ( t ) + 27 x- (t) + 27 d J—Jl x _ (t-e)+ u> 2 _ x- (t) = (■ t ) . ( 5 ) 

Here the damping term and the retarded term are derived from the nonlocal 
expressions in § , whose cutoff dependent component is absorbed with the bare 
frequency Wf, into a renormalized frequency w. The normal-mode frequency co± 
is then defined by = to 2 ± a, and the damping constant 7 by 7 = e 2 /Sttitl. 
The unit step function in Q and ([5]) clearly indicates that once the idfs come 
into interaction with the bath at t = 0, it takes some finite time £ for the 
disturbance in the field environment induced by one of the detector to reach 
the other detector, where £ is the separation between the two detectors (with 
c = h = 1). Subsequently the modified evolution of the second detector will 
prompt and send new bath disturbance back to the first one; this back and forth 
process depends on the earlier evolutionary histories of both idfs and is thus non- 
Markovian. For strong oscillator-bath coupling, measured by 2y /w±£ > 1, this 
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non-Markovian evolution can be shown to be unstable m- This behavior is in 
stark contrast to the strong coupling { n //w± > 1) Markovian dynamics, such as 
the dynamics of the CoM mode for the two idfs in the same location © where 
there is no mutual influence through the field, which tends to be overdamped. To 
avoid this instability, we will assume weak oscillator-bath coupling ( 27 /w±£ < 1 ) 
throughout. In addition, it can be shown that non-Markovianity is exponentially 
suppressed in the high temperature region (/3/£ -C 1) but merely algebraically 
in the low temperature limit (f3/£ 1). This implies that at high temperature 

(scaled by the separation), non-Markovianity is negligibly noticeable, but plays 
an increasingly important role in the system’s nonequilibrium evolution and 
entanglement dynamics at low temperatures. 

Since this retardation effect depends on the separation, if we place the two 
detectors sufficiently close to one another, we observe that the relative mode 
damps at a much slower rate than the CoM mode 033 , 


X+ + 47 X+(t) + w+ x+(t) + ■ ■ 

• = -£+(*)> 
m 

(6) 

(t) — ■*•_(*) + u 1 2 _x~{t) + ■ ■ 

• = -£-(*), 
m 

(7) 


where • • • represents the higher-order terms from the Taylor expansion of the 

2 7 

u ± - ■ 


retarded terms, and w± = or ± a =F —£-■ We see that dissipation results from 


the third-order time derivative for the relative mode, which is typically weaker 
than the counterpart for the CoM mode by the order ur±l 2 in this short sep¬ 
aration case. If this term were inadvertently excluded, then the relative mode 
would appear to be described by an undamped oscillatoiQ and the information 
about the initial state of this mode would seem to last forever, instead of being 
damped away. It would lead to a completely different dynamics of the relative 
mode. Thus these non-decaying behaviors should be more precisely understood 
as transients which last only within the time scale (yw^T 2 ) -1 . 


1 The stochastic force £_ diminishes in the short separation limit since £_=£]— £2 and 

£1, £2 are evaluated at almost the same spacetime point. 
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4. The Covariance Matrix 


The equations of motion Q, 0 enable us to compute the elements of the 
covariance matrix, defined by 


V = - Tr 
2 


p{X,X T } 


( 8 ) 


where X T = (XhPhXztPz) i n the case of a bipartite system and p is the 
corresponding density matrix. For example, the Vn element is (Xi(t)) = 
(X+W) + \ (X-(*)), where 

1 


(x±0)) = d[ ±)2 {t)( X ±{ 0 )) + —o 4 ±)2 (*)(P±(°)> 


+ — ds ds'd. 


! ^( ± ) 2 /'o'i^( ± ) 2 /V 


\s)d\ 


(a')«±006t(*')>- (9) 


Here, d^] (t) are a special set of homogeneous solutions to 0 and |5]), satisfying 
d^\0 ) = 1, 4 ±) (0) = 0, = 0, 4 ±} (0) = 1, and they are all equal to 

zero for t < 0. The noise correlation functions can be expressed in terms of the 

Gh(0 , s — s') + 
and 


Hadamard functions of the environment, (£_|_(s)£ + (s , )) = 


2 L 


G h (z 1 -z 2 ,s-s') , (£_(s)£_(s')) = 2e 2 G h (0,s-s')-Gh{z 1 -z 2 ,s-s') 
(£ + (s)£_(s')} = 0- Since the functions d^(t) in this case damp exponentially 
with time t, after the whole reduced system is fully relaxed, the terms that 
depend on the initial conditions (x±(0)), (p±(0)) will be negligible at late times. 
The other elements of the covariance matrix can be constructed likewise. 


5. Entanglement Measure 

The covariance matrix for a Gaussian continuous variable system is finite 
dimensional. This make it possible to construct the entanglement measures, 
such as negativity A f(p) and logarithmic negativity Ej/{p) [Q, GO, based on the 
smaller symplectic eigenvalues rj < of the partially-transposed covariance matrix 
V pt . Negativities are calculable, and more importantly, they offer unambiguous 
quantification of entanglement for a symmetric two-mode Gaussian state III 
Hi. They are defined by 

J\T(p) = max{0, 1 } , Etf(p) = max{0,-ln2^<} . (10) 

Zri< 
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Figure 1 : (a) Typical behavior of 77^ as a function of separation t between the two oscillators. 
An entangled state exists when 77J. < 1 / 4 ; otherwise the state is separable. There may be two 
different critical separations where the entanglement disappears, (b) Higher bath temperature 
will raise the curves upwards to make late-time entanglement harder to survive. 


Entanglement occurs when 77 < < 1/2, and the degree of entanglement can be 
described by the negativity. In the current case, the symplectic eigenvalues 77 ^ 
of ~V pt take on particularly neat forms 

1f< = min {(x+)(P-)> (x-)(p+)} , vi = max{(x+)(p?_), {x 2 J){p 2 + )} ■ (11) 

This makes interpretation of entanglement accessible. The idea is that for the 
entanglement to exist, we would like to have the uncertainties of the corre¬ 
sponding canonical variables as small as possible in order for 77 ^. < 1/4. The 
functional form of will depend on the choices of the parameters er, 7 and w± 
and its typical behavior is shown in Fig. [lj- (a). The structure of the late-time 
entanglement is much more complicated due to the interplay between different 
couplings and finite oscillator separation. 


6. Effective Description 

From the hindsights of the detailed calculations, Eqs. Q, ([5]) imply that 
at late time the normal modes can be effectively described by two uncou¬ 
pled, damped, driven oscillators with the effective oscillating frequencies W± = 


2 7 


7 


ur ± (7 =F — and the effective damping constants T± = 7 ± -- sinw±£, re- 

t UJ±t 

spectively. The non-Markovian effect is still encapsulated by the expression 
j/E. This allows us to heuristically interpret the behaviors of entanglement 

















Figure 2 : The curves of 77^, related to negativity by are plotted with respect to the oscil¬ 
lator separation for different (a) direct coupling strengths a and (b) oscillator-bath coupling 
7. Larger a and smaller 7 favor the existence of late-time entanglement. These are examples 
at T = 0. Finite temperature will raise all of these curves upwards, as shown in (b) of Fig. [T| 

based on the observations that when the temperature is not too high, the un¬ 
certainty of the normal-mode canonical variables is about the order (apart from 
the mass scale) (x±) ~ 0(W ± ), (pi j_) ~ 0(W±), and that at high temperatures 
their values are dominated by the temperature, leading to (x±) ~ 0(/3~ 1 W^ 2 ), 
(p±) ~ ©(/J -1 ). At zero temperatur^J the behaviors of rj^ are determined by 
either 0(W_/W + ) or 0(W + /W _), depending on W + ^ W_. This suggests 
that the curves (x+)(P-)> and (x~)(p+) in Fig. [j}( a ) intersect in the vicin¬ 
ity of t ~ 27 /cr, which signifies a division between the relative strengths of 
the direct inter-oscillator coupling and the field-induced non-Markovian influ¬ 
ence. When a > 2j/£, which means that the direct coupling dominates over 
the indirect field induced effect, we have r/ < ~ 0(W-/W + ). We expect that 
entanglement is possibly improved by 1) stronger inter-oscillator coupling, 2) 
weaker oscillator-bath interaction and 3) larger separation since those condi¬ 
tion may decrease the values of 77 ^ in the regime a > 2y /t. On the other 
hand, if a < 2y/£, the role of direct coupling is insignificant, and the dynam¬ 
ics is largely governed by the non-Markovian field-induced process. We have 
7 y< ~ 0(W + /W_). The same arguments suggest that entanglement will be 
enhanced by 1) weaker inter-oscillator coupling, 2) stronger oscillator-bath in- 


2 The low-temperature case applies as well because the modification in uncertainties is 
usually algebraically small and of higher orders in j 3 ~h 
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teraction and 3) shorter separation. Results involving non-Markovian processes 
are less intuitive; however, from the previous effective description at least we 
easily see how they affect the uncertainties of the normal-mode variables, which 
determine the entanglement behaviors. These qualitative descriptions based on 
detailed calculations also point out two interesting facts that direct coupling 
and the field-induced interactions between the idfs of the two detectors play a 
competing role, and that entanglement between the subsystems can possibly 
be maintained at distances much longer than previously expected, when there 
exists direct coupling between them. 


7. Critical Temperature 


At a higher temperature, the field retardation effects vanish very rapidly as 
the ratio fi/t < 1, so the dynamics of the coupled-oscillators open system be¬ 
comes simpler. In the current configuration entanglement cannot survive at very 
high temperatures because both (y;±(oo)) and (p±(oo)) are proportional to 
in the high temperature limit. Their products can be easily much greater than 
the critical value 1/4 of r)<. Thus the critical temperature for given coupling 
strengths and separation should fall in the intermediate range of the bath tem¬ 
perature. This implies neither a low- nor a high-temperature approximation can 
give an accurate prediction of the critical temperature but the high-temperature 
approximation can still offer a very reasonable upper bound for the critical tem¬ 
perature and thus is sufficient for our discussion on the critical temperature and 
its dependence on the coupling strengths. 

Here we only discuss the case a > — because stronger inter-oscillator cou¬ 
pling is a necessity to counter the thermal fluctuations/excitations from the 

bath, so as to possibly maintain the entanglement between the oscillators at 

2y 

higher bath temperatures. The other range, a < —, is less interesting since 
from the viewpoint of the effective frequency, entanglement can be sustained 
only if we require W- to be as large as possible to counter a large /3 _1 . It 


implies that 


2 7 


has to be close to 1. In such a limit the reduced system tends 
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to be unstable. Moreover, non-Markovianity is severely impeded at high tem¬ 
perature such that it becomes ineffective. In the former case, the symplectic 
eigenvalue rj < takes the form 




+ 


2 7 , P 


+ T7T 


0u\ ' 0iu\ ' 12 


1 27 , A 

~a ^- 11 7 - 1 “ ' 

0 7T Afl 


( 12 ) 


with Ap being some number much greater than /? _1 . Apparently in the 
high temperature limit increases quadratically with p^ 1 , meaning that thermal 
fluctuations dominate and introduce very large uncertainties in the canonical 
variables of the oscillators. This is particularly transparent in the large sepa¬ 
ration limit that rP. is roughly given by rfi- ~ —-- )-•••. Its value in the 

0 2 UJ+ 

high temperature limit can be brought down only if is sufficiently large. For 
fixed uj this can be achieved by increasing the inter-oscillator coupling strength, 
but only to a certain extent. The mutual influence due to separation plays a 
minor role. Larger separation only minimally alleviates the detrimental effect 
on entanglement due to thermal fluctuations. 


Suppose we extrapolate (12) to the intermediate range of bath temperature, 
and use it to identify the critical temperature by rfe = 1/4. We obtain 


Pc 


\/6 , ( VG 

-h 71 


UJ+ 
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21 A 
7 tuj± m -— 
+ A« 


+ 


(13) 


Eq. (13) gives a lower bound of 0 C , thus equivalent to the upper bound of the 
critical temperature 0~ 1 . This and the earlier qualitative analysis all consis¬ 
tently give a relation that the critical temperature should be at most about the 
order of the magnitude 0~ x = 0(u + ), that is, 


0 C GJ+ = 0(1) . 


(14) 


Eq. (131 also shows that direct coupling plays a more important role in determin¬ 


ing the critical temperature than the other factors such as the oscillator-bath 
interaction strength which shows up in 7 , the oscillator separation t and the 
field cutoff parameters. 
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8. Conclusion 


With this, we conclude that for systems well represented by two coupled 
harmonic oscillator detectors interacting weakly with a common heat bath it is 
highly unlikely that quantum entanglement can survive at high temperatures 
[Eq. @], To drive up the critical temperature one should [from Eq. ©] 
increase the direct coupling strength between the subsystems over larger sepa¬ 
rations. Both statements comply with our intuitions. This settles the question 
we raised in the beginning for nonequilibrium (Case B) systems with time- 
independent coupling in a common thermal environment. 
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